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Abstract 

A simple hadronic axion model is proposed in the framework of gauge- 
mediated supersymmetry breaking. Dynamics of Peccei-Quinn symmetry 
breaking is governed by supersymmetry breaking effects and the Peccei-Quinn 
breaking scale fpq is inversely proportional to the gravitino mass. The grav- 
itino mass range which corresponds to the axion window fpg ~ 10^ GeV 
- 10^^ GeV lies in the region predicted by gauge-mediated supersymmetry 
breaking models. The model is also shown to be cosmologically viable. 
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The Peccei-Quinn (PQ) mechanism [Q is so far the most attractive framework to solve 
the strong CP problem. The essential ingredient of this mechanism is a global U(1)pq 
symmetry which is, apart from breaking by the QCD anomaly, spontaneously broken. The 
resulting Goldstone mode (axion) acquires a small mass due to the QCD anomaly. When one 
considers the PQ mechanism in the context of supersymmetry (SUSY), there appears a non- 
compact fiat direction associated with the U(1)pq Goldstone mode. The degeneracy of the 
vacua is resolved by effects of SUSY breakdown. The properties of the PQ mechanism should 
therefore depend on details of the SUSY breakdown. In this paper, we would like to propose 
a simple model which works as the PQ mechanism in the framework of gauge-mediated 
SUSY breaking, while most of the previous study has been done in gravity-mediated SUSY 
breaking. 

We shall first describe our proposal. We consider a KSVZ axion model p[ (so called a 
hadronic axion model). A gauge-singlet Peccei-Quinn (PQ) multiplet X and also new PQ 
quarks Qp and Qp (3 and 3* in SU(3)(7) are introduced. Their U(1)pq charges are assigned 
as Q[X] = +1, Q[Qp] = -1/2 and Q[Qp] = -1/2. The superpotential of the PQ sector 
takes the following simple form: 

WpQ = XpXQpQp, (1) 

where Ap is a coupling constant. Note that no mass parameter is introduced in the super- 
potential. Eq. (|I]) gives the potential of the scalar fields in the SUSY hmit as 

V{X,QuQ2) = Xl 



2 

(2) 

where Qs and T"' denote the QCD coupling and generators. X is the PQ scalar field, and Qi 
and Q2 are PQ squarks.Q 

^Here Qi = Qp and Q2 = Qp- 



Q2Q1 



+ 1^1 



Qi + Q2 



+ 
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From the potential (|]) one finds a flat direction along the X axis with Qi = Q2 = 0. 
This comes from the fact that the superpotential (|l]) has an extended U(1)pq symmetry 
Because it is a holomorphic function, the superpotential which possesses the conventional 
U(l)pQ symmetry is also invariant under the complex form of the U(l) transformation (i.e. 
the dilatational transformation). This extended symmetry results in the fiat direction of the 
scalar potential (0) when the SUSY is an exact symmetry. However, the fiat direction is 
lifted through SUSY brealiing effects. 

In gauge- mediated SUSY break;ing theories (for a review see ref. 0]), the standard model 
(SM) gauge interaction transmits SUSY break;ing effects from a messenger sector to ordinary 
squarks and sleptons at two-loop level. At the same level the PQ squark;s Qi and Q2 also 
feel SUSY breaking through the QCD interaction, whereas the PQ scalar X is a gauge 
singlet and thus feels SUSY breaking through the Yukawa interaction ([l|) only after the PQ 
squarks receive SUSY breaking masses. Then the induced potential for X is so suppressed 
that we should not neglect the SUSY breaking effect mediated by gravity. Indeed as we 
will see shortly, the balance between the gravity-mediated effect and the gauge-mediated 
effect determines the minimum of the X field. In the following, we estimate the potential of 
the fiat direction induced by the SUSY breaking effects by considering the both mediation 
mechanisms. 

Effects of SUSY breaking communicated by the gravity are expected to induce the soft 
SUSY breaking mass to the PQ scalar field comparable to the gravitino mass ^3/2 which 
is much smaller than the electroweak scale in the gauge-mediated SUSY breaking models 
(m3/2 ~ 10~^ keV-1 GeV). Then through the gravity-mediation the fiat direction obtains 
the potential Vh as 

VH{X)^^eml/,\X\\ (3) 

where is a dimension-less parameter of the order one, which we assume to be positive. We 
have neglected higher order terms because they are suppressed by the gravitational scale. 
On the other hand, effects of SUSY breaking through the gauge-mediation mechanism 

3 



also generate a potential Vg to the PQ X field. In most models of the gauge-mediation 
SUSY is broken by some non-perturbative dynamics in a ^^idden" sector and its effects are 
fed down to a messenger sector. In the messenger sector, a gauge-singlet chiral multiplet Z 
is supposed to have a A-component vacuum expectation value (vev) (Z) and a F-component 
vev {Fz) as well. This singlet couples to messenger quark multiplets qm and which are 
3 and 3* in SU(3)c,0 in the superpotential 

Wm = Am Z qmIm (4) 

with \m a dimension-less coupling. The superpotential induces a mass of the messenger 
quark as m^,^ = M and the messenger squarks q± obtain masses squared as = 
M{M ± A). Here we define M and A by 

\m {Z) = M, (5) 
\m {Fz) = MA. (6) 

Then the SUSY is broken in the messenger sector and its effects are mediated to the ordinary 
sector through the SM gauge interaction by integrating out the heavy messenger fields. 

The PQ squarks Qi and Q2, similar to the squarks in the standard model sector, obtain 
soft SUSY breaking masses from two-loop diagrams of the messenger multiplets, gluon and 
gluino. Then the SUSY breaking effects are transmitted to the PQ X field by "one-loop" 
diagrams of Qi and Q2- 

We made an explicit calculation of the effective potential for the PQ scalar X, Vg'(^), 
along the fiat direction induced by the above gauge-mediation mechanism. Here we present 
only the final result and the details are found in Appendix.^ To estimate the minimum of 

^ In fact, we also introduce messenger leptons so that non-color SUSY particles acquire soft 
masses. 

^ After submitting the paper, we were informed that the saxion potential from the gauge-mediated 
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the potential, the derivative of the potential with respect to is sufficient rather than Vq 
itself. The result is given byQ 

dX^ (477)4 ^ [ J' 

Here we have presented the asymptotic form which is valid only for X ^ M. This is because 
for X < M, the gauge-mediation generates a sizable negative curvature for the X field in a 
standard manner, which shifts the vev of the field far away from the origin. Thus we expect 
that the minimum of X is much larger than M. This is a good news for phenomenology 
because there is some chance that the resulting PQ scale falls into the cosmological axion 
window, which is indeed the case, as we will see shortly. 

The potential along the fiat direction X generated by the two SUSY breaking effects is 
thus given by 

V,ff{X) = VH{X) + VG{X). (8) 

Using eqs. (0) and (0) we can estimate the vev of the PQ scalar field, Xc = {X) ^ M, by 
looking into the minimization condition of the potential 

which yields 

1 fas\ MA , {\lXl\ . 

The PQ scalar field thus develops a vev, and the U(1)pq symmetry is spontaneously 
broken and the PQ scale is given as /pq = Xc- We would like to emphasize that in the 



SUSY breaking effect was also calculated in ref. Q in a different manner. See section 5.2 of Q for 
more detail. Our result is in agreement with that of [0. 

^Y here should be understood as its modulus |Y|. We use this notation to avoid unnecessary 
complications. 



present axion model the PQ scale is determined by only the SUSY breaking scales without 
introducing any other mass scales. We present the numerical values of the PQ scale in Fig.l. 
If we take the parameters of the messenger sector as A = 40 TeV and M = 50 TeV the PQ 
scale is approximately given by 

f,,.X.^l^(U^), (U) 

for Ap = 1. Therefore this model provides a simple description of the PQ mechanism in 
the gauge-mediated SUSY breaking theories. Note that in the simplest structure of the 
messenger sector, our choice of A yields the right-handed slepton mass of about 90 GeV. On 
the other hand, the choice of M is somewhat arbitrary, though M must be larger than A to 
avoid a negative mass squared for the messenger squark, and thus one should regard (pi]) as 
the lower bound of the PQ scale for a given gravitino mass. One should keep this point in 
mind, though in the following we take M = 50 TeV for our representative value. 

In SUSY theories, the axion forms a chiral multiplet. In addition to the axion, which can 
be chosen as the imaginary part of the complex scalar X, the multiplet contains a saxion s, 
the real part of X, and an axino a, the fermionic component of the multiplet. Parameterizing 
X = {s/V2 + Xc)exp{ta/{V2Xc)) the mass of the saxion is easily evaluated as 

rris = ^^3/2. (12) 

Thus the saxion acquires a mass comparable to the gravitino, while the axion a gets a mass 
from the QCD effect as ~ f-nrriT,/ fpg. 

We shall next consider cosmology of this axion model to see whether or not it is cosmo- 
logically viable. The axion relic abundance from the misalignment is given as [^,0 

with 9 being the initial misalignment angle (|^^| < tt) and h the Hubble constant in units of 
100 km/sec. /Mpc. In ignorance of the misalignment angle and with theoretical uncertainty 
of order unity in the evaluation of Eq. (|TB|), we could say that the axion may be able to 
constitute the dominant component of the dark matter of the universe if 
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lO^'GeV < fpQ < 10"GeV. (14) 

The same argument implies an upper bound of the axion scale fpq < lO^^GeV. On the 
other hand, the cooling of the SN 1987A puts a lower bound of about 10^ GeV |T^. Thus 
the following range 

lO^GeV < fpQ < lO^^GeV (15) 



is regarded as the allowed region of the PQ scale. With eq. ([Tl|), this in turn identifies the 
allowed range of the gravitino mass in this scenario: 

lOOMeV > ^7713/2 > lOkeV. (16) 

Remarkably the gauge-mediated SUSY breaking naturally gives the gravitino mass in this 
range. The region of 7713/2 is somewhat shifted if we take a different value of M, but still 
the gravitino mass range fits with the gauge-mediation mechanism as far as M < 10^ GeV 
with A ~ 40 GeV fixed. On the other hand, this model does not give the possible hadronic 
axion window fpq ~ 10^ GeV. 

Let us next consider cosmology of the saxion. It dominantly decays to two axions through 
the interaction given by 

d.X^d^X = -d.sd^s + -d.ad^'a + -^—sd^ad^a + ■■■. (17) 



The decay width is estimated as 



from which we find the lifetime 

„ /lMeV\^ 1.3 X lO^sec /iMeVV 
~1.3xl0W ~ . 19) 

\ nis J ^5 V "^3/2 / 

It varies widely from 10^^-10^^ sec. for the gravitino mass range (p!6|). With this long 
lifetime, we have to see its cosmological consequences carefully. 
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To go further, we need to estimate the saxion abundance. The saxion decouples from 
the thermal bath at the temperature about W^ 



If the decoupling temperature were lower than the reheat temperature after the inflation, 
then the yield of the saxion (defined by the number density of the saxion divided by the 
entropy density) would be Ys ~ 10^^. However the gravitino problem in the GMSB suggests 
relatively low Tr. For the gravitino mass range considered here, the closure limit of the 



gravitino gives ||T^J13|1 

Tn < lO'^GeV 



VSOOGeVy VlMeVy' ^ ^ 



where rrig denotes the gluino mass. So the decoupling temperature is higher than the reheat 
temperature, and saxions are not thermalized.Q However saxions are produced by scattering 
e.g. qq — > sg in the thermal bath, with the yield |]IHP!5| 

,;.10-(^)^10-(^)-^(^). (22) 

and thus 

\ -2 



VlMeVy ^IQiiGeV; V106GeVy ^ ' 

The saxion is also produced as its coherent oscillation. An inspection shows that with the 
low reheat temperature given by (|2l| ) the saxion starts to oscillate while the infiaton obeys 
coherent oscillation (before the reheating process is completed) . We estimate the abundance 
to be 



^ When fpQ is closed to its lower bound 10^ GeV the decoupling temperature can be larger 
than the reheat temperature. Subsequent argument in the text does not apply for this case. For 
example, the abundance of the saxion simply becomes mgYs ~ 10~^ GeV (rris/lOO MeV) instead of 
(p^. We have considered cosmology of the saxion in this exceptional case separately and checked 
that this survives cosmological constraints described below. 
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where Sj is the initial amphtude of the saxion oscillation, which is expected to be of order 



In the above, we have assumed that the vacuum expectation value of the field X is 
displaced from the origin during the reheating process: otherwise thermal effects would 
generate a positive mass for X, trapping the X field at the origin until the universe cools 
down below the weak scale. In this case, the above estimate is no longer valid, and we expect 
that the relic abundance of the saxion will be much larger, which may be problematic. The 
displacement of the field X may be achieved by fluctuations of the field during the de Sitter 
expansion in the inflationary phase. In the following, we use the maximum of (^Sp and 
(^) as the saxion abundance. It remains constant as time goes on as far as no entropy 
production occurs. Note that the abundance is proportional to the reheat temperature. 

Let us next consider bounds on the abundance of the saxions and the axions produced 
by the saxion decay. Existence of exotic particles around the cosmic temperature T ~ 
1 MeV may accelerate the expansion of the universe and increase the neutron-to-proton 
number ratio at the decoupling, resulting in too much ^He abundance. Roughly speaking, 
the abundance of such exotic particles should not exceed that of one neutrino species, namely 



This is not significant. 

The energy densities of the saxions and the produced axions may alter the evolution of 
the universe at much later time, e.g. the age of the universe and the time of matter-radiation 
equality. However, if the abundance of the produced saxions does not exceed the critical 
density of the universe (divided by the present entropy density), namely 




rUsYs < 10-^GeV. 



(25) 



rrisYs < 3.6 x IQ-^h^GeV 



(26) 



then the standard cosmology is not affected. We find that the condition ( P^ ) is easily satisfied 



for a certain range of the reheat temperature. 



The next bounds we consider are those from rare decay modes of the saxion. Since 
the saxion is hght, it does not decay to gluons. On the other hand, it generally decays to 
photons. The branching ratio of the radiative decay of the saxion, B, is typically of the order 
(a/47r)^ ~ 10~^. For < 10^ sec. the constraint on the radiative decay comes from the 
photodissociation of the light elements (for a recent analysis, see ref. |TB[)- For > 10® sec, 
the non-observation of the distortion of the cosmic microwave background gives a tighter 
bound on the abundance of the produced photons: B x mgYg < 2 x 10^^^ GeV (lO^'^sec./ r^)^/^ 
forr, ~ 10^-10^° sec. and5xm,F, < 2x10"^^ GeV {lO^hec./ Ts)^/^{noh'^)-^/^ for Ts ~ 10^°- 
10^^ sec. [|17|. Comparing the above bounds with the predicted abundance (^31) or (|24D 
multiplied by the very small branching ratio B, we find that the saxion abundance easily 
survives these bounds even with relatively a high reheat temperature close to the upper 
bound allowed by the gravitino problem. On the other hand, when the saxion decays after 
the recombination (r^ > 10^^ GeV), a stringent constraint comes from the diffuse X-ray 
backgrounds. The strongest one is i? x rngYs ^ 10^^'' GeV for ~ 10^^ sec. [^. We find 
that this severe constraint is fulfilled if the reheat temperature is lower than 10 GeV so that 
the saxion abundance is highly suppressed. 

Summarizing the above arguments, we conclude that our model can be viable in view of 
the saxion cosmology. 

It is noteworthy to mention that the case where the saxion abundance exceeds (^) is 
not immediately ruled out. Rather the saxion may modify the universe's evolution in an 
interesting way. Probably the most interesting case |T9|JT^ is that the saxion energy density 
once dominates the universe before it decays to axions. Then the produced axions behave 
as an extra contribution to the radiation when we consider the cosmological developments, 
which delay the epoch of the matter-radiation equality. In this case, the saxion may play a 
role of the late decaying particle which will reconcile the standard cold dark matter (CDM) 
dominated universe for h > 0.5 with the large scale structure formation, if the following is 
fulfilled [HI 
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In particular, our model provides a concrete realization of ref. , where a hypothetically 
light saxion is proposed to play this role. 

Finally we will look into the axino, the fermionic superpartner of the axion. Radiative 
corrections with the messenger fields give an axino mass. In our case, the contribution 
appears at the three-loop level, and it is ~ 777,3^2/^ ^3/2- The axino mass is also induced 
by supergravity effects and in general it becomes of the order of the gravitino mass 



An inspection shows that, however, in this simple model, the leading contribution of the 
order 7773/2 is cancelled, and the non-vanishing contribution is ~ (/pQ/M^)7r73/2 ^ ^3/2- 
Thus the axino is much lighter than the gravitino. The axinos are produced by scattering 
after the inflationary epoch and we expect that the yield of the axinos is more or less the 
same as that of the saxions (|2^) . In our case, since the axino is much lighter than the saxion, 
the axino energy density is always smaller than the saxion's one and thus it is harmless. 
Another possible worry is that the gravitino decay to an axino and an axion is kinematically 
allowed. However, with the mass range of the gravitino considered here, the lifetime is much 
longer than the age of the universe, hence the gravitino is essentially stable. Thus the light 
axino appeared in this model does not alter the argument above and it is cosmologically 
harmless. 

To conclude, we proposed the hadronic axion model in the gauge-mediated SUSY break- 
ing theories. The model gives a simple description on the PQ symmetry breaking mechanism 
which is solely governed by the SUSY breaking physics. It should be noted that the PQ 
scale is inversely scaled to the gravitino mass and the resulting PQ scale naturally falls into 
the axion window, for the gravitino mass range favored by the gauge-mediation. Moreover 
the axion becomes a dark matter candidate for 7773/2 ~ 10 keV-1 MeV. Next we investigated 
cosmological implications of the other particles in the axion supermultiplet, namely the sax- 
ion and the axino, and showed that our model is cosmologically acceptable with sufficiently 
low reheat temperature which is suggested by the gravitino problem. Finally we pointed out 

11 



that the saxion in this model may play a role of the late-decaying particle, causing the delay 
of the matter-radiation equality, and consequently reconciling the standard CDM with the 
structure formations of the universe. Details of this issue will be discussed elsewhere ||23 . 
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APPENDIX: EFFECTIVE POTENTIAL FOR PQ SCALAR FIELD X 



Here we derive the potential for the PQ scalar field X induced by the SUSY breaking 
effects through gauge-mediation mechanism. We calculate an effective potential for Xc, the 
classical value of the X field. As the ordinary squarks, the PQ squarks receives the SUSY 
breaking effects from the messenger sector by the QCD interaction at two-loop level. Then 
the effects are transmitted to X by the Yukawa interaction (|I]) at three-loop order. The 
Feynmann diagrams for the vacuum energy which should be estimated are shown in Figs. 
2, 4 and 5. In our calculation we expand it in terms of the MA and estimate the leading 
term of (9((MA)^), i.e., the leading term with F-vev of the messenger multiplet. From the 
Yukawa interaction (P the PQ quark and squark acquire a mass Mq = XpX^. 

First, we estimate the vacuum energy induced by the D-term scalar potential which 
diagram is shown in Fig. 2. This diagram gives 

d^k f d'^p f d^q 



iV, 



D 



(27r)4y (27r)4y (27r)4 
1 

[k"^ — m~ ] [{k + pY — m~ ] 



1 



[q^-Ml][{q + pf-Ml]\ 



(Al) 



where the messenger squark masses are =M(M ± A). Then we expand the RHS in 
terms of (MA) as 

, r d^k f d^p f d^q , /"i 
^^"='^^7 (2^7 (2^7 (2^ io ^"7o 



(27 

1 



+ 



g2 + y(l-y)p2_M2 

3(l-2a;)2(MA)2 



+ ■ ■ 



(A2) 



\k^ + x{l - x)p^ - M^Y ' [k^ + x{l - x)p^ - M^f 
where the (9 (MA) term vanishes on making the integration on x. Since what we are inter- 
ested in are the SUSY breaking effects at the first non- vanishing order, we pick up only the 
0{{MAy) term. Thus the effective potential from the D-term interaction is estimated as 

d'^k f d'^p f d^q 



iVd = gt 



■1 rl 

(27r)4y (27r)4y {2n)^ Jo '^'^ Jo '^^ 
12(1 -2x)2(MA)2 



g2 + y{l - y)pi - + x(l - X)p2 - M^f 



(A3) 
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Next, we turn to the vacuum energy with the gluon fields. We first calculate the one- loop 
correction to the vacuum polarization from the messenger multiplets, (see Fig. 3.) and then 
estimate the vacuum energy shown in Fig. 4. Up to the 0{{MA)^), the result is given by 



^^^ = ^^7 (2^7 (2^7 (2^70 ^^7o 

36(1 -2x)2(MA)2 



g2 + y{l - y)p2 - [P + x(l - x)p'^ - M'^f 

Finally, the vacuum energy induced by the gluino (see Fig. 5.) is estimated as 

d'^k f d^p f d^q 



(A4) 



-9t 



I dx dy 
Jo 



(27r)4 7 (27r)4 7 {2tt)^ 

192(1 -y)(l -a;)3(MA)2 



(A5) 



g2 + y{l - y)p2 - [P + x(l - x)p'^ - M'^f 

Therefore summing up eqs. (|A3|) , ([AiD and ([A5|), we obtain the effective potential Vq at 
the order (MA)^ as 

d^k r d^p f d^q /■! ^ /■! 



VG{MX)=tgt{MKr 



{2'kYJ {2tiYJ (27r)4 
48x(l - x) 



I dx dy 
Jo 



g2 + y{l - y)p'i - M^] + x(l - - M"^]" 

and by making a Wick rotation it leads to 

Vg{MI) = ^ (MA)' jdkjdpjdq £ dx £ dy 

x(l — x)kpq 



(A6) 



(A7) 



q + y{l-y)p + M^ [k + x{l - x)p + M'^]'^ 

Instead of the effective potential itself, it is more convenient to calculate its derivative. We 
find 



a;(l — x)kpq 



q + y{l-y)p + Ml [k + x{l - x)p + M^f 
Now the momentum integrations are straight forward. 



(A8) 
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dx / dy 



X 



{MAY 
(47r)6 M2 Jo 

■\ — x{l — x) + x(l — x) ln[z^x(l — x)] 



— x{\ — x\ 



where 



M, 



Q 



M2y{i-y) • perform the x integration analytically. (See ref. 

here we give an approximate expression valid for z"^ ^ 4 

\2 



aVb(M2) 8gi {MAY 



(4vr)6 Ml Jo 

Sgt {MAY 
{ArcY Ml 



dy In 



M'r 



M^y{l-y) 



M2 



The PQ squark mass is expressed by the vev of X and we obtain the final result 



dXt 



9t (MA)'j^2 



(477 



X2 



M2 
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FIG. 1. Values of the PQ breaking scale fpQ. Parameters for the messenger sector are listed 
in the figure. 
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FIG. 2. The vacuum diagram induced by the D-term scalar potential. q± are the messenger 
squarks and Qi,2 are the PQ squarks. 




FIG. 3. The one-loop diagrams of the vacuum polarization induced by the messenger multiplets. 
q± are the messenger squarks and denotes the messenger quark. 
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FIG. 4. The vacuum diagrams with gluon. Qi^2 are the PQ squarks and *q denotes the PQ 



quark. A blob represents the one-loop vacuum polarization shown in Fig.3. 




FIG. 5. The vacuum diagram with gluino g. 
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